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Abstract 

The spinless Falicov-Kimball model on the simple cubic lattice is analyzed 
in the Hubbard-I and dynamical mean field (DMFT) approximations. The 
Matsubara and real frequency itinerant electron Green's functions, the evolu- 
tion of the system with doping, and the range of phase separation are found 
in two approximations. At large values of the on-site Coulomb repulsion 
both approximations give similar results. The phase separation can be also 
favorable for a more general model, where heavy electrons have a finite band- 
width. This indicates that the phase separation phenomenon is an inherent 
feature of the systems described by the Hubbard-like models with wide and 
narrow bands. 



1. Introduction 

Strongly correlated electron systems, the systems where the interaction 
between electrons is larger than their kinetic energy, are intensively stud- 
ied nowadays. Heavy fermion compounds [if, high T c cuprates P(, magnetic 
oxides such as manganites j^, cobaltites, and many others, belong to this 
class of materials. In theoretical study, the Hubbard model and its modifica- 
tions such as t — J model and others are widely used to explain the complex 
behavior of strongly correlated materials. The main problem arising in the 
analysis of the Hubbard-like Hamiltonians is the absence of reliable analytical 
methods at large electron-electron interaction. The numerical studies such 
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as Monte-Carlo simulations and exact diagonalization on small clusters have 
large computational costs. On the other hand, there exists a huge number 
of different approximate methods not requiring time-consuming calculations. 
However, each of these methods has a restricted range of applicability, and it 
is not always clear whether the results obtained by this or that approximation 
reveal an actual physics or just reflect the features of initial assumptions. 

The simple example of the models of correlated electrons is the Falicov- 
Kimball model. It was proposed in seminal paper of L. M. Falicov and 
J. C. Kimball, Q, to describe the metal to insulator transition in rare-earth 
and transition-metal compounds. It includes both itinerant and localized 
electrons with on-site interaction between them. This model exhibits many 
interesting features such as non Fermi-liquid behavior ^ for any value of 
coupling constant U, charge-density wave phase 0,0,0], and phase separa- 
tion js), lil, [lO , 11 ] . The Falicov-Kimball model has been studied both numeri- 
cally with Monte-Carlo simulations and exact diagonalization SB [3, and 
by using different approximations jsl, Hi, 15]. Several exact results have been 
obtained also, such as the existence of the phase separation at large coupling 
constant 10 



The phase separation is a quite ubiquitous phenomenon in the strongly 
correlated electron systems. It manifests itself in different situations (both 
in experiment and theory), where the itinerancy of charge carriers competes 
with their tendency to the localization. However, revealing the phase separa- 
tion in a given system or model is not a trivial task, which usually demands 
a lot of numerical and analytical work. The first most reasonable step is the 
using of the analytical mean-field type approximations, which do not require 
laborious computations. If we reveal the tendency to the phase separation 
in the framework of these approximations, a problem arises, whether this 
result is an artefact of the approximation, or it captures the main physics 
and remains qualitatively unchanged in a more refined (but much more cum- 
bersome) approximations such as dynamical mean-field theory (DMFT) or 
in numerical simulations. 

To address this problem, in this paper I use as a testing ground the 
Falicov-Kimball model for which the DMFT results are quite transparent 
and compare DMFT with the Hubbard-I approximation. More precisely, I 
analyse the Hubbard-I and DMFT approximations of the spinless Falicov- 
Kimball model in three dimensions in the strong U limit. The Hubbard-I ap- 
proximation |l6[ first proposed for the Hubbard model can be applied to the 



Falikov-Kimball model as well 111 . This is a simple mean-filed-like scheme, 
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which allows to obtain analytical results. Although this approximation fails 
for U < W, where W is the bandwidth of itinerant electrons, it can be con- 
sidered as a good approximation at large U. A more advanced approach, dy- 
namical mean- field theory |l7| , has been applied first to the Falicov-Kimball 
model in jfj| , and has been considered later in many theoretical papers (for a 
review, see, e.g., [Isj)- This approximation takes into account the local corre- 
lations more rigourously, but at the same time it requires more computation 
costs. Here I show that for strong coupling, U ^> W, both approximations 
give similar results for the band structure of itinerant electrons, the evolution 
of the system with doping, and the region of the phase separation. 

After this, I analyze the possibility of the phase separation for a more 
general model, where 'localized' electrons have a finite bandwidth Wf. This 
model can be considered as a two band Hubbard model for spinless fermions. 
The itineracy of 'localized' electrons makes DMFT calculations much more 
cumbersome in comparison to that for the Falicov-Kimball model. But again, 
both approximations give similar results for the range of the phase separation. 
It is shown that the phase separation for this model exists when Wf < Wf, 
where W* ~ although the Hubbard-I and DMFT approximations 

give slightly different values for the WJ. This correlates well with the result 



of [19|, where the phase separation in the spin- 1/2 two-band Hubbard model 
is shown in the Hubbard-I approximation. So, the results presented in this 
paper indicate that the phase separation phenomenon is the inherent feature 
of multiband Hubbard models with wide and narrow bands, and not an 
artefact of this or that approximation scheme. 

2. Falicov-Kimball model in the Hubbard-I approximation 

The Hamiltonian of the Falicov-Kimball (FK) model is Q 

H = _J2(tv4dj+H.c)-e f Y,n{ + uY,n*n{ 

ij i i 

+ nf), (1) 

i 

where nf = d\di, n{ = fjfi, and d\ (dj) and fj (fa) are creation (annihilation) 
spinless fermion operators for the band and localized electrons, respectively. 
The first term in (JTJ) describes the kinetic energy of the band electrons, the 
second term controls the position of the localized level with respect to the 
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center of the conduction band, while the third term corresponds to the on- 
site interaction between the band and localized electrons. In further analysis 
I will consider attractive interaction with U > 0. 

The Hubbard-I approximation is based on the decoupling in the sys- 
tem of equations for the electrons Green's functions. The use of Matsub- 
ara Green's functions is more convenient to study the system at finite tem- 
peratures. It also allows simply comparing the Hubbard-I and DMFT ap- 
proximations, since DMFT deals with Matsubara one-particle Green's func- 
tions. The equation of motion for one particle <i-electron Green's function 
G d (T-T',t-j) = -(fd t (r)d](T')) is: 

(—Iff + A G ^ 1 -j) + Y1 1 - i) ~ 

UF d (T,i-j) = 6 tJ 5(T), (2) 

where F d (r — r', % — j) = — (Tn{ (r)di(r)dt(r')) is the two-particle Green's 
function. The equation of motion for F d is the following: 

d \ 
-- + H-U) F d (r,i-j) 

- 5>,(T n{ (rH(r)4(0)> = n% 8(r) , (3) 
i 

where n* is the averaged number of the localized electrons. In the Hubbard-I 
approximation, the following decoupling in the second term in ([3]) is made [16( : 

r, % — j). Under this substitution, the closed 
system of equations for the G d and F d is obtained. As a result, one ob- 
tains the following expression for the one-particle Green's function (in the 
frequency- momentum representation) : 

n (■ i \ iu n + fi-U{l -n f ) 

G d {iu n ,k ) = - ■ — — jT-yr , 4 

[lU n + fl - E + (k)\ [lUJn +11- E-{k)\ 

where u n = (2n + 1)ttT are the Matsubara frequencies, 



£±(k) .£(^ ± ^£M_^) 2 + £ ,„/ £(k) , (5) 

and e(k) = —^2jtije k( - ri ~ r ^ is the spectrum of free (i-electrons. The band 
structure of interacted (i-electrons consists of two Hubbard subbands with the 
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spectra ^(k). The filling of lower (upper) subband is equal to 1 — rz/ (w). 
The gap between Hubbard subbands exists for any value of U. This is the 
main weakness of the Hubbard-I approximation: it can not explain the metal 
to insulator transition existed in the Hubbard-like models. The one-particle 
Green's function of the /-electrons Gf(r — r', i — j) = —(Tfi^fUr')) is 
obtained in the same manner. The result for Gf in the Hubbard-I approxi- 
mation is 

_ v 1 — n d n d 

Gf{iu n , k) = , + . (6) 

iu n + n + e f iU) n + fi + ef — U 

The /-electrons Green's function remains local, but the localized level splits 
into two sublevels corresponding to sites occupied by or free of ci-electrons. 
The filling of the lower (upper) level is equal to 1 — n d (n d ). 

The knowledge of one-particle Green's functions allows to evaluate the 
thermodynamics quantities of the system. The average number of <i-electrons 
is given by the relation n d = TJ2 n J dk/VBze luJn0 Gd(iUn,k), where Vbz is 
the volume of the first Brillouin zone. Taking the summation over Matsubara 
frequencies, after straightforward algebra, one obtains 



n d 



where £"™ ax and E± m are the maximum and minimum values of the functions 
£±(k), 

Mm = ^ + 1 , fi = ~, (8) 
and po(E) is the density of states of free (i-electrons 

£L6(E-eQ£)). (9) 

VBZ 

This function is non-zero only in the region \E\ < W/2, where W is the 
bandwidth of free <i-electrons. Taking into account that po(E) = outside 
this region, one can extend the integration over E in ([7]) to all real axis, and 
omit the summation over two Hubbard subbands. 
The average number of /-electrons is 



n 



f = (1 - n d )n F (-e f ) + n d n F {U - e f ) . (10) 
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When the total number n of electrons is fixed, the chemical potential is deter- 
mined from the relation n = n d (fi) + n* (//). This equation, and equations (!7j) 
and (TlO|) set the system of equation for the self-consistent determination of 
fi, n d , and rJ . For the purposes of the present study, it is more convenient to 
consider the system at fixed value of rJ . In this case, the chemical potential 
fi is found from ffTUl) . At small temperatures, /i becomes fixed at the point 
fi = —€j when n < 1, or at // = — e/ + U when n > 1. 

The density of states for (interacted) (^-electrons in the Hubbard-I ap- 
proximation is 



v £- (1 -n-QC/, 
The free energy of the system F can be written in the form 



dEp(E)$(E) + 

(1 _ n d M-e f ) + n d ^(U - £f ) , (12) 
where 

$(E) = (i n F (E) - T In [l + e ^~ E) ] . (13) 

The one-particle Green's functions and the mean values n d and w as func- 
tions of n will be compared in the next sections with those calculated using 
DMFT approximation. 



3. Falicov-Kimball model in DMFT approximation 

In the dynamical mean field approach, the lattice problem is replaced 
by the effective single-site problem. The interaction of electrons in a given 
site, say i — o, with other sites in the lattice is described by the effective 
dynamical field A(r) which has to be determined self-consistently. For the 
Falicov-Kimball model, the single-site Hamiltonian is 

H = -e f n f + Un d n f - p(n d + n f ), (14) 

where n d = Sd, n* = pf, and d\ /' (d, f) are creation (annihilation) 
operators for d- and /-electrons at site i = o. The partition function Z for 
this single-site model can be represented as a functional integral of exp(— S&) 
with the effective action S c s = L dr{d)dd/dr+pdf /8t+Hq) + Jq dr J Q dr' 
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g^(t)A(t — r')d(r'). The partition function can be also written in the inter- 
action representation with the non-local in time interaction operator Ui Qt (r) 
= J c£t / g^(t)A(t — T')d{r'). Both representations are, of course, equivalent. 
In general, for the Hubbard-like models, the partition function can be cal- 
culated only numerically since H is not quadratic in electronic operators. 
Fortunately, for the Falicov-Kimball model, it is possible to obtain an ana- 
lytical formula for Z. This is related to the fact that n* commutes both with 
H and U int . 

In the interaction representation the partition function has a form 



Z = Tr [e-^a] 



where 



and 



a = f exp 
d(r) 



- dr dT'd){r)\{T -r')d(r') 
. Jo Jo 



(15) 
(16) 



THo de- rHo 



d}(r) = e THo Se- THo . (17) 

Since the operator n* commutes both with d{r) and H , the trace over / 
operators is trivial and gives 



Z = Tr 



_|_ e £(M+e/)Tr 



(18) 



where traces are taken over d operators only. Expanding each term in per- 
turbation expansion series, taking the traces by applying the Wick theorem, 
and re-summing the series again, one obtains: 



Z=[l + e^]Z x {y) + e^ +e f\l + e^~ u) ]Z x ^ - U), 



(19) 



where 



OO s 



A, 



iu n + E 



(20) 



and A n = f Q dre JWnT A(r) is the Fourier transform of the function A(r). 
The local <i-electrons Green's function is 



1 5Z 

Z 5\ n 



n 



f 



n J 



iu n + [i — \ r 
1 



+ 



(21) 



iu n + fJL — Xn-U iu n + fM- A n - £„ 
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where E n = T,(ico n ) is the self energy of the effective single-site model, and 
w is the averaged number of /-electrons 



n 



f 



T5Z _ 
ZST f ~ 

I _|_ e -^(M+ e /). 



(22) 



l + e 



n -1 



Z A (/i - [/) 1 + ePb-u) 
The average number of (^-electrons is 

- Wnr ■ " 1 1 



n 
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(23) 



The c?-electrons Green's function in the lattice can be written in the form 

1 



G d (iu n , k) 



/i-e(k) - E(iu;„, k) 



(24) 



where E(iu; n , k) is the self energy of the lattice model. In the DMFT ap- 
proximation, it is assumed that the self energy £(zu; n ,k) does not depend 
on the quasimomentum and equals to E(zo;„,k) = E n , with £„ from ( 1211) . 
This is uncontrollable scheme which becomes exact, nevertheless, in infinite 
dimensions 17 . From (EH), one obtains 



D(ioj r , 



dk 



V, 



G d (iu n , k) = D (iu n + jj, - E„) 



BZ 



where 



W/2 



D{E) 



dE 



, Po(E') 
E — E' 



(25) 



(26) 



-W/2 



Equations ( 1211) and (1251) allow to calculate self-consistently both A„ and 
E n . From ( 12T|) and (j25|) one obtains 



A, 
E, 



iu n + n - E n - — 

D {icu n + fi- E n J 



n 



iuj n + ji — \ n — (1 — nf)U 



(27) 
(28) 



where n/ is given by (1221) . This system of equations is solved by iterations. 
In the iteration scheme it is more convenient to work either at fixed number 
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of localized electrons w or at fixed total number of electrons n, instead of 
fixed chemical potential \i (especially at low temperatures). For given A n , 
E n , and /i, the new values of A n are calculated from ff27|) . the new value of 
/i is found from f )22|) if n* is fixed or from the relation n = n d (fi) + n^(/j,) if 
n is fixed, and, finally, the new values of E n are evaluated from ff28l) . For 
the relative accuracy 10~ 4 , this iteration procedure converges in about 10 
iterations. 



4. Comparison of Hubbard-I and DMFT approximations 

Consider first the itinerant electrons Green's function. All results below 
are obtained for the simple cubic lattice with the spectrum of free (^-electrons 
of the form e(k) = — H / (cos/c :r + cos/c ?; + cos/c z )/6. In the DMFT approxima- 
tion the (i-electrons Green's function is given by f]2~4"]) . where E(ia; n , k) = E n 
is calculated from the system of equations (122]) , ff27|) , and ([28]) • The typical 
curves of real and imaginary parts of E n vs u n are shown in figured] for small 
(U = 0.1 WO and large (U = 5W) values of U. From (T57D and Q2B) it is 
possible to obtain an analytical expression for E n in the narrow band limit, 
when W — > 0. Expanding the function D(E) in the Taylor series up to the 
second order, from the system (1271) and (128]) one obtains: 



A„ 



24 



\ — rJ nf 
+ 



iuj n + (j, iuj n + jM — U 



+0(|r), (29) 



iw n + /i-(l-n/)?7' 1 1 



Substituting f]30]) into (]24|) . after simple algebra, one arrives to @ for the 
Green's function in the Hubbard-I approximation. The self energy cor- 
responding to the Hubbard-I approximation is also shown in figure [H It is 
seen that at small U the asymptotic (130]) has qualitatively different behavior 
in the region uj n < W than E n calculated from (127]) and (128]) . For large U 
both curves practically coincide even for small values of u n . 

Thus, both approximations give similar results for large U. Note, how- 
ever, that limits U ^> W and W — >■ are different, and even at U — > oo there 
will be a discrepancy between E n and E^ at small frequencies u n W. One 
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Figure 1: The dependencies of Re£ n and Im£ n vs oj n calculated for U/W — 0.1 (upper 
panel) and U/W — 5 (lower panel). Other parameters are: rJ = 0.4, ef/W = 0.1, 
T/W = 0.01. Solid (dashed) curves correspond to the DMFT (Hubbard-I) approximation. 



can say, that the Hubbard-I approximation corresponds to the DMFT one 
in the narrow band limit W — > 0. 

The use of Matsubara Green's functions is convenient for calculation of 
the thermodynamics quantities. But if we are interesting in the spectral prop- 
erties, we need working at real frequencies. There are two possible ways for 
transition from Matsubara to real frequencies. One of them is the analytical 
continuation of calculated Matsubara Green's function to real axis. It can be 
done by different methods 17, 2(3, 21]. In this paper I use Pade approximants 
method 21J in evaluating the real-frequency Green's function. Another way 
I propose in this paper is the use of analytical continuation of the system 
of equations ff271) and ff28l) which is obtained by replacing iu n —> u + iO. 
The infinitesimal positive imaginary part is added to calculate the retarded 
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self-energy. The system of equations (127)1 and ( 1251) becomes 



D R [lo + 11 - S(w)) 

S M = nfU TT n77> ( 32 ) 

where D R (E) = D(E + iO), and for real E, D R (E) = D(E) - iTrp (E). The 
values of \x or/and rv in (13TT) and (|32|) are found by solving the iteration 
scheme at Matsubara frequencies as described above. Excluding the A(a;) 
from the system (I5Tj) and (13"2"j) . one arrives to the equation for the self-energy: 

- u] = i ( , 

E(u)-nfU D R {u + fi-E{uj))' 1 ' 

This equation is solved numerically. Such method of finding the retarded 
self-energy is more reliable than the Pade approximants method of analytical 
continuation of S n , especially at high frequencies, where Pade approximants 
method depends on the cut-off frequency and can fail in some cases. But at 
frequencies \u\ < W both methods give the same results. 

It is possible to obtain an analytical expression for S(cj) in two different 
limits: for small U, and for small number of /-electrons, n' < 1. In the 
former case, expanding A and E in (!3l!) and (1321) in series of U, one obtains 
in the second order: 

S(w) = n f U + n f (l -n f )U 2 D R (u + n) , (34) 
ImE(u) = -7m f (l-n f )U 2 p (co + fx). 

This result coincides with E calculated by perturbation theory in the second 
order of U. Note a non- Fermi liquid behavior of the self-energy 0, Et| : 
Im E(0) 7^ even for T — > 0. This is due to the infinite mass of /-electrons. If 
we allow of /-electrons to have a finite bandwidth Wf, then the perturbation 
theory gives the standard Fermi-liquid result for the imaginary part of the 
self-energy at frequencies \co\ >C Wf 

ImE(w) oc -(w 2 + T 2 ). (35) 

At w < 1, the expansion of A and E in the series of n* gives in the 
leading order: 
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Figure 2: The frequency dependence of the spectral function (|37|l at different values of e 
and for three values of U: U/W = 0.1, U/W = 0.3, and U/W = 1 for panels from top 
to bottom. Other model parameters are: n' = 0.4, tf/W — 0.1, and T/W — 0.01. In 
the lowest panel, circles connected by dashed lines mark the position of the 5-peaks of the 
spectral function calculated in the Hubbard-I approximation. 



This is a well known result for the self-energy of the electrons scattered on 



localized impurities in the limit of small concentration of impurities [22 
Thus, in both limits, the DMFT approximation gives results confirming by 
controllable approximations such as perturbation theory. 

Consider now the spectral properties of the conduction electrons. The 
spectral function of ci-electrons is A(u,k) = — ImG^a;, k)/7r. Since in both 
approximations all momentum dependence contains in the non-interacted 
spectrum e(k), it is more convenient to consider the function 

A(u>, e) = -- Im P ° (g) - , . (37) 
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Figure 3: The density of states in DMFT (solid curves) and Hubbard-I (dashed curve) 
approximation calculated at different values of U. Other model parameters are: w = 0.5, 
e f = 0, T/W = 0.01. 

The curves of A(u,e) vs u at different values of e are shown in figure [2] for 
three different values of the interaction constant U. At small U, there is only 
one peak in the spectral function for any value of e. In an intermediate case, 
the two-peaked structure of the spectral function appears at some range of e. 
Finally, at large U, there exists two separate Hubbard bands of <i-electrons. 
The self-energy in the Hubbard-I approximation has no imaginary part, the 
spectral function consists of two 5-peaks located at frequencies u = E±(e)—p, 
where E± are given by (J5j) with e(k) = e. The position of these peaks are 
shown by circles in the lowest panel of figure [2] corresponding to U — W. It 
is seen from this figure, that both approximations give very similar results 
for the band structure of (^-electrons in the large U limit (with the exception 
of finite lifetime of quasiparticles in DMFT). 

The (i-electrons density of states p(E) is obtained by integration of A(E — 
p,e) over e. Using (|33|) . from (l37j) one obtains: 

, . 1 T / 1 - nf n f \ 

p(E) = - Im — r + — . 38 

HK ! tt \E(E - p) - U E(E - p) J y 1 

The density of states (138]) calculated for the symmetric case, e/ = 0, n f = 0.5, 
and at different values of U are shown in figure |3j The density of states ( ITT]) 
in the Hubbard-I approximation is shown for U = W. In this case, the lower 
Hubbard subband turns out to be filled, n d = 0.5. It is seen, that in the 
DMFT approach the metal-insulator transition occurs at U c = 0A5W (in 
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Figure 4: The dependencies of n d and n* on n. Solid (dashed) curves correspond to 
the DMFT (Hubbard-I) approximation. The parameters of the model are: U/W — 1.5, 
ef/W = 0.1, T/W = 0.01. The vertical arrows show the position of doping n when the 
chemical potential equals to the bare and renormalized energy of the localized level (see 
the text below). 

the symmetric case, e/ = 0, and for simple cubic lattice). Below, I will study 
the case U > U c , where both approximations give qualitatively the same 
results for the density of states. 



5. Phase Separation 

In the previous sections, the homogeneous state of the Falicov-Kimball 
model was studied. But it is well known that for wide range of the model 



parameters an inhomogeneous state is preferable at low temperatures [18 
It can be either charge density wave, or phase separation. For the symmetric 
case (ef — 0) and at half- filling, the checkerboard phase is realized jsl, 0, S], 
while away of half-filling the phase separation prevails jsl, 0, 0, [ll| • 

In this paper I study non-symmetric case (Ef ^ 0) and restrict myself 
by consideration of the phase separation only. To do this, I consider the 
evolution of the system with doping and examine the stability of the homo- 
geneous state by analyzing the system's free energy. The typical dependencies 
of numbers of d- and /-electrons, n d and n*, on doping n are shown in fig- 
ure HI Both approximations give qualitatively the same behavior of n d and 
nf vs n. At small doping, n < 1, only <i-electrons exist in the system since 
the chemical potential lies far below the localized level. With the increase 
of n the /-electrons appear. But the filling of correlated band depends on 
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Figure 5: The free energy of the system vs doping. Solid (dashed) curve corresponds to the 
DMFT (Hubbard-I) approximation. Dotted curve is the free energy of free <i-electrons. 
In the region where d- and /-electrons coexist the free energy has negative curvature 
indicating the instability of the homogeneous state. Maxwell construction to both curves 
is shown by dotted-dashed lines. The parameters of the model correspond to figure [4j 



the number rJ of localized electrons: both the bandwith and the filling of 
the lower Hubbard band decrease when n* increases. Thus, in the further 
increase of n (and tj/), the number of itinerant electrons starts to decrease. 
Finally, at n ~> 1, all localized levels become occupied, and (^-electrons in the 
upper Hubbard band appear. 

The quantitative disagreement between dependencies n d (n), n/(n) is re- 
lated mainly to the different behavior of rv on the chemical potential /i in two 
approximations. For the Hubbard-I the dependence of w on \i is governed 
by (ITU]) . In the case U ^> W and /i < U, it can be approximated as 

, 1 - n d 
n f « -. - , . 39 

Thus, at small temperatures, /-electrons appear when chemical potential 
reaches the value \l — — £/■ In the DMFT approximation, the dependence of 
n> on fi is given by ( I22p . When 1, one obtains from the system (1271) 

and (128]) : A n ~ iu n + \i — 1/D(iu n + /i). At T — > 0, one can replace the 
Matsubara frequency summation by integration over u, and as a result, one 
obtains: 

U ~ e -/?[e} (/*)+/*] + I ' ^ 40 ^ 
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Figure 6: The phase diagram of the Falicov-Kimball model at U/W — 3 and Sf/W = 0.1. 
Inside the closed region the phase separation exists. Solid (dashed) curve is calculated in 
the DMFT (Hubbard-I) approximation. 



where 



and 



't(fi) =e f - /„(//, n) + f (/j, - U, fj) 



f (E,/i) = Ee(E) + 



In 



7T 



E - [i + l/D(iu + /j) 



E 



(41) 



(42) 



Thus, in the DMFT approximation, the interaction between electrons renor- 
malizes the energy of the localized level [5(. The /-electrons appear when the 
chemical potential reaches the value defined by the equation /x + £*(//) = 0. 
Since the renormalized level has larger energy, the DMFT approximation 
gives the larger value of doping n when /-electrons appear (and, consequently, 
the doping when n d starts to decrease, see figure HJ). With this exception, 
both approximations give the same behavior of n d and n* vs n. 

The free energy in the Hubbard-I approximation is given by ffT2"j) . In the 
DMFT approximation, the partition function (fl9|) defines the grand potential 
of the effective single-site model: Qi oc = —Tin Z. The grand potential Q and 
the free energy F = Q + [in of the lattice model can be found following the 
approach described in 0. As a result, one obtains: 



F = Q ioc + fin + T^2 de p {e) In 



lUJr, 



-fi-e-Ti n 



(43) 
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The dependence of the system's free energy on doping calculated in two 
approximations is shown in figure EJ The model parameters correspond to 
figure HI The free energies in two approximations coincide in the regions of 
n where < 1 and at n ~ 1 (where n' i < 1). In the region of doping 
where itinerant and localized electrons coexist, the free energy has negative 
curvature indicating the instability of the homogeneous state. In this region 
the system separates into the phases with fi^ k and w ~ 1. With the 
increase of temperature, the doping range n\ < n < ri2 where phase sep- 
aration exists becomes narrow and disappears at some critical temperature 
T c . This result is summarized in the phase diagram presented in figure |6j 
The Hubbard-I approximation gives the value of n\ lower than the DMFT 
one. This is related to the fact that in the DMFT approximation /-electrons 
appear at larger doping due to renormalization of the energy of the localized 
level as described above. In addition, the DMFT gives the lower critical 
temperature. 

6. Extension to the Hubbard model 

Commutativity of the Falicov-Kimball model's Hamiltonian with the num- 
ber of /-electrons strongly simplifies the DMFT scheme for this model; it is 
possible to obtain an analytical expression for the single-site partition func- 
tion Z. For other models of strongly correlated electrons, such as Hubbard 
model, the DMFT approach requires much more CPU time. The partition 
function Z and the local Green's functions can be calculated either numeri- 
cally or by using some approximations which work only in the special ranges 
of model parameters. 

In this section I consider a more general model where /-electrons have a 
finite bandwidth. The model Hamiltonian is the sum of the Falicov-Kimball 
model's Hamiltonian (TIJ, and the term describing the kinetic energy of the 
/-electrons: 

#/ = -^£(/M + i? - c -)' ( 44 ) 

{ij) 

where Wf is the bandwidth of free /-electrons. The k-dependence of the 
/-electrons spectrum is assumed to be the same as for (i-electrons. Such a 
model can be considered as a two-band Hubbard model for spinless fermions. 

In the Hubbard-I approximation, the expression for Gd{ioo n , k) remains 
the same, but the /-electrons Greens's function, Gf, and the dispersion (for 
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Figure 7: The occupation numbers, n d and w , vs doping at finite bandwidth of /-electrons. 
Solid (dashed) curves correspond to the DMFT (Hubbard-I) approximation. The param- 
eters of the model are: W f /W = 0.2, U/W = 2.5, e f /W = 0.1, and T/W = 0.01. In the 
inset the system's free energy is shown. Solid (dashed) curve corresponds to the DMFT 
(Hubbard-I) approximation. The Maxwell constructions for both curves are shown by 
dot-dashed lines. 



interacted electrons) take the form of (j4]) and (jSJ), respectively, with the 
replacement wJ — > n d , \i —)■// + ef and e(k) — >■ e(k)Wf/W. 

The changes in the DMFT scheme are more significant. Now one needs to 
introduce two dynamical mean-fields, \d(iu> n ) and \f{iu n ). These fields, and 
the local Greens's functions, Dd(ioJ n ) and Df(iu n ) [together with self-energies 
Yidiiujn) and T,f(iu:„)], are calculated numerically by the 'exact diagonaliza- 
tion' method 17|, |23|]. In the iteration scheme, the total number of electrons 
is fixed, and for each iteration the chemical potential is recalculated from 
given \ d j{iu n ) and E d>f (iuj n ). 

The detailed comparison of two approximations for such this model is 
beyond the scope of the present paper. Here, I restrict myself by the analysis 
of the stability of the system toward the phase separation. The mean values 
n* and n d , and the system's free energy F as functions of doping can be 
calculated in the way similar to that described in Sections [2] and |3j The 
typical curves are presented in figure [71 The value of the /-electrons band 
is taken to be Wf = 0.2W, and the energy e/(= 0.1W) characterizing the 
relative position of the bands corresponds to figure HI Both approximations 
give qualitatively the same behavior of occupation numbers and the free 
energy, very similar to that presented in figures H] and |5j The quantitative 
difference can be explained again by the renormalization of the energy e/ 
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appeared in the DMFT approximation and lost in the Hubbard-I one. In both 
approximations the system turns out to be unstable toward separation into 
the phases with rJ ~ and n d « 0. The region of phase separation depends 
on the value of Wf. it narrows when Wf increases and disappears at some 
critical value Wf (which depends on other model parameters). The Hubbard- 
I approximation gives Wf/W ~ 0.3 for the model parameters corresponding 
to figure while the DMFT approximation gives smaller value, WJ /W ~ 
0.25. 

7. Conclusions 

In conclusion, I compared the results giving by the Hubbard-I and DMFT 
approximations for the spinless Falicov-Kimball model in the simple cubic 
lattice. The itinerant Green's function, the evolution of the system with 
doping, and the possibility of the phase separation were considered. It was 
shown that both approximations coincide in the narrow band limit, W — > 0. 
For the finite W and the strong coupling, U ^> W, both approximations do 
not contradict to each other giving qualitatively similar results in the spectral 
properties of itinerant electrons and the behavior of the system with doping. 
Both approximations predict the instability of the homogeneous state toward 
the phase separation in the wide range of model parameters. It was also 
shown that the system remains unstable toward the phase separation even 
for finite bandwidth Wf of heavy electrons. This indicates that the phase 
separation phenomenon is the inherent feature of Hubbard-like models with 
wide and narrow bands, and not an artefact of this or that approximation 
scheme. 
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